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We introduce variants of the Ma-Dasgupta renormalization-group (RG) approach for 
random quantum spin chains, in which the energy-scale is reduced by decimation built 
on either perturbative or non-perturbative principles. In one non-perturbative version of 
the method, we require the exact invariance of the lowest gaps, while in a second class 
of perturbative Ma-Dasgupta techniques, different decimation rules are utilized. For the 
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matrix renormalization-group (DMRG) calculations and previous studies. 
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1. Introduction 

In recent years intensive research work has ex plored some curious properties of 
antiferromagnetic Heisenberg chains .^J^MIEIB] ^pj-^g features of the pure systems 
are well-known after Haldane's seminal wor kin half-integer 5" and integer S chain 
show distinct behavior; gapless spectra, quasi-long-range order and gapped spec- 
tra, hidden topological order characterize ther n, r espectively. In S = 1/2 chain. 
Ma and Dasgupta introduced a RG techniquej^'^ with which it was possible to 
point out that any small amount of disorder triggers the random-singlet phase 
(RSP) where spins in arbitrarily large distance are coupled in singlets. In S" = 1 
chain,''' 8 9 10 11 12 13 14 15 Q^^y sufficiently large disorder is able to change the 
pure system properties. 

The S = \ random antiferromagnetic Heisenberg chain (R AHC) was studied nu- 
merically with density-matrix renormalization-group (DMRG) ^^^ l ^^ ^ ^^ l ^'^ l ^^ l ^^ and 
quantum Monte Carlo (QMC) technique, a considerable part of the find- 
ings was contributed by extended versionJ^tSEHOl Qf the Ma-Dasgupta RG method 
that is based on a perturbative decimation step, which fails for higher-S RAHC 
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Jl J2=i2 J3 J' 

Fig. 1. The four-spin system with the largest coupUng parameter Q replaced by a two-spin system 
with coupling constant given in Eq. (|3ll, in a usual SDRG decimation, or such J' which involves an 
identical first gap as that of the original four-spin subchains, in a GRG decimation. 

for weak disorder, but asymptotically correct for strong disorder, therefore, it is 
known as st rong-disorder RG (SDRG). The SDRG method has been applied in 
chains ,^2ZI28| ladders, ■^^ in two-, and three-dimensional systems,-^ for a review by 
Igloi and Monthus, see Ref. [511 

In this paper, a non -perturbative RG technique is presented, which is inspired by 
the SDRG methocPSl and in which the lowest gaps play the central role, therefore it 
is termed as gap RG (GRG). The original SDRG is investigated: a simple argument 
is presented and justified, by using local and global aspects of disorder, to explain 

fool 

the success of the SDRG, which holds true in higher-S" chains.- ■ 

The structure of the paper is the following: in Section 2, the model and the 
failure of the SDRG decimation are described. The results of the GRG and the 
disorder-induced phases are described in Section 3. Different SDRG methods, the 
interpretation and justification of their correctness, are detailed in Section 4. The 
results are summarized and compared to earlier findings, in Section 5. 

2. The random antiferromagnetic Heisenberg chain and the 
strong-disorder renormalization 

We investigate the S — 1 RAHC with Hamiltonian 

L-l 

H = ^ JiSi ■ Si+i. (1) 

2=1 

The Ji > couplings follow a random distribution, the power-law distribution 

ps{J) ^S-\j-^+^/^ forO<J<l (2) 

where the S parameter tunes the strength of disorder as 5^ — var[lnJ]; S — 0, 
S ~ 1.0, and S ^ 00 denote pure, uniform, and i nfin ite randomness, respectiv ely. 
This system has a relevance in solid state physics,'S3l in quantum comput ations 
and serves as a testing ground for developing novel methods and theories .^^^^^ 

In a given RAHC, the SDRG method eliminates iteratively the largest coupling 
with its four neighboring spins and replaces it by two spins with a new coupling, 
see Fig. [U determined by perturbation theory for general S as 

^, _2 S{S+l)JiJ3 

^ - 3 n ■ ^"^^ 

The method works in S = 1/2 RAHC where the new coupling J' — JiJ-i/2Vl 
is smaller than any of the eliminated ones. For higher-S' RAHC the (2/3)5(5 -I- 
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Fig. 2. Integrated probability distribution of the first gaps at (5 = 0.5 (a) and at 5 = 1.0 (b) 
disorder. In terms of the slopes of the low-gap region, the agreement is obvious, although the GRG 
data (empty symbols) are shifted toward the low-gap region. DMRG data (filled symbols) are 
available up to L = 32, the GRG data are plotted also up to this length. The estimated slopes are 
indicated in the legends (i.e., {DMRG slopes} and [GRG slopes]) and plotted with dashed lines. 



1) > 1, the new couphng can be larger than the decimated one: the method fails 
for wea k disorde r . In order to circumvent this problem, novel RG schemes were 
initiatedMini 

that use a set of decimation steps and deal with a mixture chain 
of S = 1/2 and S = 1 spins. These RG schemes predict a transition fro m the 
pure system behavior to the RSP, which was tested with DMRGII and QMC. 1221141 
Saguia et al^^ use a different strategy to tackle the problem, they keep the original 
SDRG decimation, if it works correctly, and replace it with a two-step perturbative 
decimation, if it fails. The method yields three phases: gapped/gapless Haldane 
phase (or nonsingular/singular Griffiths phases), and RSP. 



3. Gap renormalization-group method 

We briefly describe a non-perturbative RG method that has a similar iterative 
strategy as the SDRG: it eliminates the largest couplings, but the decimation rule 
is different, the new coupling is determined with the condition that the first gap 
is invariant, see Fig. [1] The system properties can be drawn from the first gap 
distribution and by analyzing the RG flows. The GRG works in higher-5 RAHC 
and in S = 1/2 chain in its original form.*^^ Preliminary results of GRG were 
published in Refl32l In the following, more detailed investigations are presented. 



3.1. The disorder-induced phases and finite-size corrections 

We have investigated the 5 = 1 RAHC by using the GRG technique and 
DMRG^"'^ ■^^ method in parallel. We eliminated iteratively L — A spins from 
the original L by using the GRG method, the new couplings were always smaller 
than the eliminated ones, and the final four-spin system is diagonalized by the 
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Fig. 3. The I/2:' exponent. Empty symbols (dashed lines) denote the GRG data (fitted straight 
lines) . Full symbols (solid lines) denote the SDRG data. Opaque symbols (dot-dashed lines) denote 
the XX chain data. Error bars, if not indicated, are smaller than the symbols. 

Lanczos method. On the other hand, genuine first gaps are calculated by using 
DMRG method. The basis number in DMRG calculations was so large (maximum 
180) that it did not influence the probability distribution. 

The comparison was made for different strength of disorder, always a very good 
agreement was found between the original gap integrated probability distribution 
and the GRG data in the sense that the slopes of the low-gap tail from GRG and 
DMRG are the same, see in Fig. [S] In fact, this is the first time that RG results of 
5 = 1 RAHC are faced to genuine probability distributions. 

The dynamical exponent z can be determined via the scaling of the first gap 
probability densities (and integrated probability densities) in the low-gap region: 

P(lnA) - A^/^', (4) 

where z' is identical to the dynamical exponent z, if z' > 1, otherwise z — \. 
From the low-gap tail of the integrated probability distributions the slopes (1/z') 
in the log-log plot are extracted for different sizes and b parameter, Fig. [3l In this 
parameter region, the following relation is conjectured and utilized in the figure 

l/z'((5, 00) helps to identify the phase diagram of the 5 = 1 RAHC. At very weak 
disorder, up to = 0.4(15), z'(5, 00) is smaller than 1 indicating the nonsingular 
Griffiths phase. At intermediate disorder 0.4 < 8 < 1.0, z'((5, 00) has still finite 
value, singular Griffiths phase, while above be = 1.00(15) we have found z' = cxd in 
the large-i limit indicating the appearance of RSP. 
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Fig. 4. The phase diagram. 5 = denotes the pure system, <5i the phase boundary between the 
nonsingular and singular Griffiths phases, <5c the phase boundary between Griffiths phase and RSP. 




Fig. 5. Analysis of the GRG flow at <5 = 1.0. The tit exponent from the density of active spins is 
identified from the slope as 1/2. Inset: The a parameter can be identified as a = — l/lnil. 



We have performed SDRG also on XX chains and analyzed the gap distributions 
in the same size region. The XX chain is chosen because its decimation rule,l^with 
1 prefactor in Eq.®, is the closest to the decimation rule of 5 = 1 RAHC, with 
prefactor 4/3. The slopes, extracted from the XX chain at 6 — 0.3 and S = 1.0 and 
plotted in Fig.[niscale with effective exponent —1/4 and predict disappearing l/z', 
as it is expected,!^ and thus justify the assumption that these system sizes provide 
physically conclusive results in the large-L limit. 



3.2. Gap renormalization flow in the RSP 

In RSP during the RG treatment of a certain random chain, the ratio of the re- 
maining spins that are not decimated out follows the scaling 

la#] 

11 

Z " |lnr!|i/^- 

This relation describes the connection between the characteristic length scale of the 
system L and the energy scale fl via the ijj = 1/2 exponent, a characteristic relation 
of RSP. The GRG method is applied on very large system sizes, 20 million spins, 
results are plotted in FigO The slope of the double logarithmic plot reveals the 
presence of RSP in the system at (5 = 1.0 disorder. 
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Fig. 6. (Color online) The gap distributions generated by the original SDRG (empty symbols) and 
the modified SDRG (filled symbols) for <5 = 0.3, 1.0,3.0 disorder in (a), (b), and (c), respectively. 



The following aim is to show out that the scaling of the coupling distribution 
under GRG procedure follows the expected RSP behaviorP 

with a = — 1/ \nfl. Actually, the above type of scaling is found at 6 — 1.0, see the 
inset in Fig|5] where a, determined from the slope of the coupling distribution, is 
depicted as a function of — l/lnfi. Both quantities, the coupling distribution and 
the fraction of undecimated spins follow the scaling that is characteristic of the 
RSP. 



4. Results of the SDRG 

In this section, different aspects of SDRG method are analyzed and discussed. 
Firstly, proper interpretation of the strong-disorder limit is presented that explains 
the reliability of SDRG method. Secondly, the scaling relations and the correspond- 
ing finite-size corrections of SDRG data are discussed. Finally, an unusual represen- 
tation of the SDRG flows is established that has a very high accuracy. 

4.1. The low-gap region of the probability distributions 

One has to clearly distinguish two different kinds of appearance of disorder: global 
and local appearances. On the one hand, the global disorder describes the statistical 
ensemble of couplings in the chains and it can be parameterized, in this work, with 
the S parameter of the power-law distribution. On the other hand, the local disorder 
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or subchain disorder characterizes only one particular four-spin subchain and it can 
be parameterized, for instance, with the value max(Ji, J^)/Q,. 

The SDRG can be correct in the proper strong-disorder limit: if the subchain 
disorder is strong everywhere in a chain that is already a sufficient condition. In 
certain chains the strong subchain or local disorder can be present generally at 
arbitrary location in the chain with some small probability that itself depends on 
the global disorder that describes the statistical ensemble of the couplings. These 
probably rare chains provide the smallest gaps to the gap distributions since in 
these chains repeatedly very small new couplings are replaced in the consecutive 
decimations resulting the smallest final couplings and gaps. Consequently, the low- 
gap regime of the results can be correct, in RG sense, at any given strength of global 
disorder, as they come via the correct decimation steps from those rare chains in 
which the subchain disorder is strong everywhere. This interpretation of the lo w-gap 
region is in accordance with the well-established scenario of rare region effects. I^^ESl 



4.2. Comparison of the original and the modified SDRG 

In order to make this feature more transparent, we present a systematic investigation 
and comparison of SDRG and modified SDRG results in the low-gap region of 
the probability distributions. The original SDRG method is straightforward, the 
modified SDRG applies the same decimation rule except in those cases in which 
the new coupling J' is larger than at least one of Ji and J3. In those decimation 
steps cJ', with c < l/[{2/3)S{S + 1))] = 3/4, is used as a new coupling ensuring 
always smaller new coupling than the decimated ones. The two schemes are almost 
identical, the only difference is factor c if the perturbative decimation fails. 

This modified version of SDRG has a similar scenario as the RG method applied 
in Ref. However, the modification is very simple here, just a rescaling of the 
coupling given by the original rule, whereas a two-step degenerate perturbative 
decimation isis applied in Ref ll2l Nevertheless, the heuristic argument used here to 
explain the success of SDRG applies also for the modified version SDRG in Ref. [121 

These two RG schemes, the original SDRG and the modified one have been 
applied to generate the gap distribution for several system sizes and strength of 
disorder. The two types of gap distributions are plotted in Fig. \6\ the slopes of the 
low-gap regions are identical for fixed sizes and strengths of global disorder. The 
plotted results were generated with c = 2/3 in the modified decimation steps, but 
several other values provided identical slopes in the low-gap region. 

The slope of the low-gap region is apparently independent of those decimation 
steps where the original decimation rule fails in perturbative sense. The modified 
RG steps, by generating smaller coupling than the original one, shift the corre- 
sponding low-gap distributions towards the lower gaps, but the structure of this 
regime remains the same as in the original case. A similar relation of the GRG and 
DMRG results is observed in Fig.O For longer chains that involve more decimation 
steps, the shift between the two gap distributions is more enhanced, see Fig. |6l For 
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Fig. 7. (a) Gap distribution generated by the original SD RG method at (5 = 0.5. (b) Scahng 
collapse assuming Griffiths phase with scaling relation in Eq.l(9]l and z = 5.0. 

sufficiently strong global disorder the modified decimation steps practically do not 
make change on the structure of the gap distributions, at least for smaller chain 
sizes, see Fig. [6l^b) and (c) for 6 = 1.0 and S = 3.0 disorder, respectively. 

4.3. Disorder-induced phase diagram 

The original SDRG results were carefully analyzed and the estimates of 1/z' were 
extracted for different system sizes and plotted in Fig[3]as a function of i"^/"*. These 
numerical values are in agreement with the results of GRG method. The estimates 
for the phase boundaries 6i = 0.29(6), 6c = 1-0(1) are in rather good agreement 
with other studies. 

One has to notice, the finite-size corrections are smaller in the Griffiths phase 
than in the RSP, in terms of the dynamical exponent z, and are weakening deeper 
into Griffiths phase. The SDRG results are more reliable than the GRG results. 
Firstly, more random chains are averaged for SDRG, at least 80 million, and only 2 
million for the GRG. Secondly, larger systems are analyzed for the SDRG (L = 1024 
and only L = 128 for GRG). Even larger system sizes were tested, however for those 
sizes the convergence in the low-gap region was not satisfactory although intensive 
computations were performed on a Beowulf Linux cluster in order to have accurate 
gap distributions from GRG, SDRG, DMRG methods. 

4.4. Scaling of the SDRG data 

In the following, the scaling of the SDRG results in the S = 1 RAHC are tested by 
the expected scaling laws ■ in the Griffiths phase and in the RSP that are shown out 
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exactly, for instance, in the random transverse- field Ising spin chain, in the 5 = 1/2 
RAHC or in the random XX chain. In the Griffiths phase, the energy scale J7 is 
continously reducing during the RG, and the couplings follow the distribution 

Due to the relation of energy scale and the length of the system { ft ^ and by 
assuming that the RG process leaves the z exponent unaltered, one can conclude 
easily Eq.(I3|) and the scaling law 

FL(lnA) = P(ln(Ai-)). (9) 

In the RSP, where z is infinity, one can set formally z — — In A ^ L"^ , Eq. © 
transforms into Eq. ([7]) , and the corresponding scaling relation is 

FL(ln A) - i-^P(ln(A)/i'^)) (10) 

with -0 = 1/2. These s cali ng laws were success fully utilized, for instance, in the 
random XX spin chains,!^ in the S* = 1/2 RAHCJS] 

and also in higher-dimensional 
systemSjISSl but were never tested in the S — 1 RAHC due to the known reasons. 

On the one hand, these scaling laws assume that the investigated sizes are free 
of finite-size corrections. If this condition is not accomplished, and it is not accom- 
plished in the S = \ RAHC, one can expect the appearance of finite-size corrections 
in addition to the above detailed leading behavior. On the other hand, Eq.® as- 
sumes the monotoniicc decrease of the energy scale VL that is surely not rigorously 
true during the original SDRG in the S — \ RAHC, therefore Eq.® is surely not 
accurately true. However, Eq.® can be true in a looser sense, i.e., its consequences, 
Eq.®, Eq.® , Eq.® and Ea. (fTO| . can be true since the extraction of these relations 
does not assume the rigorous accurateness of Eq.®. The probability distribution 
in the low-gap region follows the expected behavior, Eq.®, as shown in the earlier 
subsection. The scaling behavior in the Griffiths phase, Eq.®, and in the RSP, 
Eq. (fTO| . are investigated in the following paragraphs. 

The SDRG gap distributions are plotted in Fig[7|for 6 = 0.5 disorder. One can 
notice a systematic broadening in the distributions and systematic change in the 
slopes of the low-gap region with increasing systems sizes, which can be considered 
as a result of finite-size corrections. The presence of finite-size corrections is evident 
also from Fig. [3l where the estimated slopes are plotted. Due to these finite-size 
corrections the expected scaling collapse is not perfect, see Figllljb). It is to be 
emphasized that the values of the z exponent in the scaling collapse are taken 
from the large-L estimates in Figl3] and they are not optimized for the scaling 
collapse itself. In order to have even better scaling collapse by using the same large- 
L limit estimates of z or to see a satisfying fit with an analytically estimated gap 
distributions,^'^ one should study drastically larger systems sizes. 

The scaling collapse is tested also in the RSP. The gap distributions of the 
original SDRG are clearly broadening with increasing system sizes. Fig. ®a), the 



March 3, 2009 8:41 WSPC/INSTRUCTION FILE lajko 



10 Peter Lajko 





-5 

sr -10 

c 

-15 



-20 

-200 -100 -8 -4 -15 -10 -5 

x=lnA x=lnA/L"' x=lnA/L"^ 

Fig. 8. (a) Gap distribution generated by the original SDRG method at 5 = 1.0. (b) The scaling 
collapse assuming RSP and the scaling relation in Eg. (1101 with ip = 1/2, and with ip = 1/3 (c). 

changes in the slopes are more pronounced than in Fig. [71 indicating broadening 
without hmit and the presence of RSP. The scahng collapse is not perfect for the 
investigated sizes, although there is apparently a crossover phenomenon present: for 
larger sizes the collapse is almost perfect, see Fig|51[b). 

In addition, the scaling collapse is examined with and exponent '0 = 1/3 as this 
value of exponent was theoretically predicted in Ref. [7l [8] and numerically found 
in Ref. [32] at the critical point, see Fig[8ljc). It can be argued that the scahng 
collapse is far better with f/) = 1/2 than with = 1/3. For ip = 1/2, the collapse 
is almost perfect in the large-gap region and there is a systematic improvement for 
larger sizes in the low-gap reg ion . Very recently, '0 = 1/3 was found with scaling 

lOOl 

investigation of DMRG results^ at the same disorder for relatively small system 
sizes. The results presented here in Fig[8] imply that the 1/3 value of exponent is 
only an effective value due to a crossover phenomenon. Since there is a uncertainty 
regarding the value of 6c, the scaling collapse of RSP was investigated at several 
different strength of disorder in the S ~ 0.7 — 1.5 region. In this region, the 1/2 
exponent gives always better scaling collapse than the 1/3 exponent. 

We think that the presence of crossover phenomena, discussed here in detail, 
associated with changing finite-size corrections is the explanation why the 5* = 1 
RAHC represents a very challenging problem often leading to controversial results 
regarding the location of phase boundaries, see Section 7.5. in Ref. [31] or Table II 
in Ref. [32] This phenomenon is reported in detail in the 5 = 1/2 RAHC in Ref. l42l 
1431 1441 but also known in other quantum spin chains. 

To summarize the scaling collapse investigation, scaling collapse considerations 
should be carefully inspected as there are affected by finite-size effect, while the 
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Fig. 9. (Color online) Symbols denote the largest coupling applying the original SDRG method. 
Thick (red) staggered lines denote the same when the modified SDRG is used. Renormalization 
flows are depicted for 6 = 0.1,0.2,0.3,0.4,0.7,0.8,0.9,1.0 in Fig. (a)-(h), respectively. The -1.3 
and —2.2 exponents in the label of the horizontal axis are determined in a fitting procedure. 

finite-size corrections allow a rather accurate estimates of different phases in terms 
of the slopes of the distributions in the low-gap regions. 

4.5. Analysis of SDRG flows and the i/j exponent 

SDRG and modified SDRG results were confronted also to relations in Eq.® and 
^ with the same conclusion as it is plotted in FiglS] at 5 = 1.0 disorder, but 
unfortunately also in a region below this strength of disorder due to considerable 
errors. This indicates that probably other representations of RG flow are more 
proper to describe the phase diagram and the appearance of the RSP. It is found 
that the best quantity to be investigated is ln(| In ri|)/ Inp as a function of In p. 
From Eq. ([6]) , it is evident that assuming the presence of RSP this product should 
tend to —0.5 in the small-p limit. 

The results of this investigation are plotted in Figl5] for both the original and 
modified SDRG cases. The data are independent of the original length of the chain, 
which is 1 million in each plotted case, they are practically free of finite-size correc- 
tions and thus allow a very precise estimate on ijj. 

For weak disorder, very large fluctuations dominate the data, FiglQl which grad- 
ually decrease for stronger disorder. Around S — 0.2 — 0.4 an envelope curve appears 
below the fluctuating data points that allows a straight line flt for small-p estimate. 
With increasing disorder the small-p estimates gradually reach the —0.5 value, at 
5 = Sc, and remains there for stronger disorder. This representation provides a very 
accurate prediction: RSP appears at Sc = 1.00(3) disorder. Notice that the small-p 
estimates gradually decrease to —0.5 and then with further increase of the strength 
of disorder they remain there, indicating clearly that the phase transition point is 
at 6c = 1.00(3) where tpc = 0.50(1) but tp = 0.5 also deep in the RSP. 
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5. Summary of the results and comparison to earlier findings 

In this paper, a non-perturbative GRG scheme and different variants of the per- 
turbative SDRG method give a coherent picture about the disorder-induced phase 
diagram of 5 = 1 RAHC. 

The GRG method provides an accountable phase diagram: for weak disorder, 
the system is in a Griffiths phase with a disorder-dependent dynamical exponent; 
for stronger disorder the systems reaches the RSP, a.t 6 — Sc, where the dynamical 
exponent is infinity. The Griffiths phase can be further divided, at S — 6i, into 
nonsingular and singular phases. GRG provides reasonable estimates on the phase 
boundaries by applying finite-size scaling analysis, see FiglHl Si — 0.40(15) and 
6c — 1.00(15). From a theoretical point of view, the GRG scheme is well-established; 
from a numerical point of view, the data from GRG technique is well-justified. 

The original and modified SDRG properly describe, in contrast to the earlier 
expectation, the features of the system and provide a high accuracy for the phase 
boundaries in the large-i limit: Si — 0.29(6) and 5c — 1.01(3), in quantitative agree- 
ment with the latest works. The exponent ■0 is shown out with different approaches 
and consistently found to be 0.5 in the RSP and in the critical point. 

The findings of this work are at least in qualitative agreement with all of the 
earlier studies. This agreement can be considered, beyond the direct DMRG test at 
smaller sizes, as indirect test of the SDRG method in the asymptotic regime. 

Good quantitative agreement can be observed with the latest numerical investi- 
gations. The 5i and Sc phase boundaries, determined in this work, are very close to 
the values in Ref. [121 [HI HSl [32l Hida's early but hidden estimate of Sc, see Fig. 3 
in Ref. [Til is also close to 1.01(3) that is the best estimate for Sc in this work. The 
phase boundary Si — 0.29(6), deep in the Griffiths phase, well coincides with inde- 
pendent estimates: Si = 0.23 (the relation of the power-law and box distribution is 
detailed in Ref. [42|) in Ref. [Hand = 0.45 in Ref. |3l 

There is a group of studies that successfully describe the phases in the S — 1 
RAHC by decomposing the S — 1 spins into 5 = 1/2 spins in an SDRG scheme. 
The success of these well-established methods is not questioned by the success of 
the original SDRG method; however, with these earlier works only a qualitative 
agreement is found, because these works contribute only qualitative description 
of two phases RSP and a Griffiths 

phasfPI 

or they predict rather distinct phase 
boundaries: Sc ~ 0.77 is found in Ref. [7l M a nd S c ~ 1-5 in Ref. [TOl 

One of the earliest extended SDRG studj^^ predicts distinct ip exponent in the 
critical point (V'c = 1/3) and in the RSP {ijj = 0.5). In the present work, a novel 
representation of the SDRG allows an accurate numerical determination of Sc and 
ip exponent in the RSP as well as at the critical point, in both cases ip = 0.50(3) 
is observed, which is supported also by scaling collapse of the gap distributions at 
and around Sc = 1.01(3) . ^ = 0-5 w as found at 5 = 1.0 also in a series of recent 
QMC and SDRG studies.H^I^IIIIlSj 

Thus, the properties of the RSP at the critical 
point and for different strength of disorder seem to be coherent as it is indicated in 
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Ref. [26]by means of QMC. 

The success of the SDRG approach, the fact that the perturbatively incorrect 
decimation steps are not influencing the low-gap region of the gap distribution, is 
an important contribution to our understanding of SDRG procedure. Furthermore, 
these accurate results are not in contradiction with the conclusion of Ref. [45l i.e., 
that the SDRG is not suitable to describe physical quantities like the free energy. In 
this work, the study of the GRG and SDRG flows provide the physical conclusion 
from the asymptotical behavior of the low-gap region within a usual RG framework 
without explicite relevance on physical quantities like free energy. Actually, this is 
obvious in SDRG context, see Ref. [3ll[3aor|46l|47l|48l|49l[50l 

It is very likely that some of these new findings can be used to study the features 
of higher-S* RAHC and also other systems.'^ It would be instructive to test in what 
extent the phase diagram for 5* = 1 RAHG is reproducible by large scale QMC tech- 
nique in terms of physical quantities that are inaccessible by SDRG. Moreove r, sim- 
ilar phases in the S = 1/2 ferromagnetic-antiferromagnetic alternating chain j^l^^ 
ladders, and whether the extended versions of the SDRG, involving S = 1/2 
spins, can be numerically refined and can provide quantitatively consistent results 
that remain subject of further studies. 
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